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We study a class of metric-variation f(R) models that accelerates the expansion without a cos- 
mological constant and satisfies both cosmological and solar-system tests in the small-field limit of 
the parameter space. Solar-system tests alone place only weak bounds on these models, since the 
additional scalar degree of freedom is locked to the high-curvature general-relativistic prediction 
across more than 25 orders of magnitude in density, out through the solar corona. This agree- 
ment requires that the galactic halo be of sufficient extent to maintain the galaxy at high curvature 
in the presence of the low-curvature cosmological background. If the galactic halo and local en- 
vironment in f(R) models do not have substantially deeper potentials than expected in ACDM, 
then cosmological field amplitudes |/h| > 10 -6 will cause the galactic interior to evolve to low 
curvature during the acceleration epoch. Viability of large-deviation models therefore rests on the 
structure and evolution of the galactic halo, requiring cosmological simulations of f(R) models, and 
not directly on solar-system tests. Even small deviations that conservatively satisfy both galactic 
and solar-system constraints can still be tested by future, percent-level measurements of the linear 
power spectrum, while they remain undetectable to cosmological-distance measures. Although we 
illustrate these effects in a specific class of models, the requirements on f(R) are phrased in a nearly 
model-independent manner. 
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I. INTRODUCTION 

Cosmic acceleration, in principle, can arise not from 
dark energy — a new, exotic form of matter — but rather 
from a modification of gravity that appears on large 
scales. The addition of a non-linear function of the Ricci 
scalar R to the Einstein- Hilbert action has been demon- 
strated to cause acceleration for a wide variety of f (R ) 
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What is less clear in the literature is whether any pro- 
posed metric- variation f(R) modification can simultane- 
ously satisfy stringent solar-system bounds on deviations 
from general relativity as well as accelerate the expansion 
at late times [H H M, M, M, S H M, M MaM, Izi • 
Chiba [75[ showed that the fundamental difficulty is that 
f(R) gravity introduces a scalar degree of freedom with 
the same coupling to matter as gravity that, at the back- 
ground cosmological density, is extremely light. This 
light degree of freedom produces a long-range fifth force 
or, equivalently, a dissociation of the curvature of the 
space-time from the local density. As a result, the metric 
around the Sun is predicted to be different than is im- 
plied by observations. This problem has been explicitly 
proven to exist for a wide variety of f{R) models, if the 
Sun is placed into a background of cosmological density 
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If high density could be re-associated with high curva- 
ture this difficulty would disappear. The scalar degree of 
freedom would become massive in the high-density solar 
vicinity and hidden from solar-system tests by the so- 
called chameleon mechanism [78|, [7^, H(| ■ This requires 
a form for f(R) where the mass squared of the scalar is 
large and positive at high curvature (rJlj . Such a condi- 
tion is also required for agreement with high-redshift cos- 
mological tests from the cosmic microwave background 
(CMB) [UH^I. It should be considered as a necessary 
condition for a successful f(R) model; it is violated in the 
original inverse-curvature model and many other gener- 
alizations (e.g. Hl|). 

Faulkner et al. [8(| analyzed a class of models where 
solar-system tests of gravity could be evaded, but only 
at the price of reintroducing the cosmological constant 
as a constant piece of f(R) that drives the cosmic accel- 
eration but is unrelated to local modifications of gravity. 
Moreover, for these models to satisfy local constraints, all 
aspects of the cosmology are essentially indistinguishable 
from general relativity with a cosmological constant. 

In this Paper, we introduce a class of f(R) models that 
do not contain a cosmological constant and yet are ex- 
plicitly designed to satisfy cosmological and solar-system 
constraints in certain limits of parameter space. We use 
these models to ask under what circumstances it is pos- 
sible to significantly modify cosmological predictions and 
yet evade all local tests of gravity. 

We begin in fJTT] by introducing the model class, its ef- 
fect on the background expansion history and the growth 
of structure. We show that cosmological tests of the 
growth of structure can, in principle, provide extremely 
precise tests of f(R) gravity that rival local constraints 
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and complement them in a very different range in cur- 
vature. We then analyze local tests of gravity in S jllll 
and show that solar-system tests alone are fairly easy to 
evade, provided gravity behaves similarly to general rela- 
tivity in the galaxy. However, if cosmological deviations 
from general relativity are required to be large, the lat- 
ter condition is satisfied only with extreme and testable 
changes to the galactic halo. We discuss these results in 

m 

II. f(R) COSMOLOGY 

In this section, we discuss the cosmological impact of 
f(R) models of the acceleration. We begin in §11 Al bv in- 
troducing a class of models that accelerate the expansion 
without a true cosmological constant but nonetheless in- 
cludes the phenomenology of ACDM as a limiting case. 
We then describe the background equations of motion 
( mi B[) and their representation as an equation for the 
scalar degree of freedom f ^II C[) . Finally, we calculate 
the expansion history ( §11 D| and linear power spectrum 
( §11 E|) in our class of f(R) models. 
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FIG. 1: Functional form of f(R) for n = 1,4, with normaliza- 
tion parameters ci, C2 given by |/ro| = 0.01 and a matching to 
ACDM densities (see S]II D|) . These functions transition from 
zero to a constant as R exceeds m 2 . The sharpness of the tran- 
sition increases with n and its position increases with |/ijo|- 
During cosmological expansion, the background only reaches 
R/m 2 ~ 40 for |/ao| <C 1 and so the functional dependence 
for smaller R/m 2 has no impact on the phenomenology. 



A. Model 

We consider a modification to the Einstein-Hubert ac- 
tion of the form [84 



S = 




(1) 



where R is the Ricci scalar, which we will refer to as the 
curvature, k 2 = 8irG, and C m is the matter Lagrangian. 
Note that a constant / is simply a cosmological constant. 
We work in the Jordan frame throughout this paper. 

We choose the functional form of f(R) to satisfy cer- 
tain observationally desirable properties. Firstly, the cos- 
mology should mimic ACDM in the high-redshift regime 
where it is well-tested by the CMB. Secondly, it should 
accelerate the expansion at low redshift with an expan- 
sion history that is close to ACDM, but without a true 
cosmological constant. Thirdly, there should be sufficient 
degrees of freedom in the parametrization to encompass 
as broad a range of low-redshift phenomena as is cur- 
rently observationally acceptable. Finally, for the pur- 
poses of constraining small deviations from general rela- 
tivity with cosmological and solar-system tests, it should 
include the phenomenology of ACDM as a limiting case. 

These requirements suggest that we take 



lim f(R) = const. . 



lim/(tf) 



0, 



(2) 



which can be satisfied by a general class of broken power 
law models 



f(R) 



j Cl (R/m 2 ) n 
c 2 {R/m 2 ) n + 1 ' 



(3) 



with n > 0, and for convenience we take the mass scale 



K 2 Pq 



(8315Mpc) 



0.13 



(4) 



where po = /5(lna = 0) is the average density today. c\ 
and C2 are dimensionless parameters. It is useful to note 
that 



2 

^4 = !- 22 8 x 10 30 . 
777. \ lg cm J 



0.13 



(5) 



The sign of f(R) is chosen so that its second derivative 



fun 



d 2 f(R) 
dR 2 



> 



(6) 



for R 3> to 2 , to ensure that, at high density, the solu- 
tion is stable at high-curvature [6l|. This condition also 
implies that cosmological tests at high redshift remain 
the same as in general relativity (GR). For example, the 
physical matter density Sl m h 2 inferred from the CMB us- 
ing GR remains valid for the f(R) models. As such, m is 
a better choice of scale than H since it does not vary for 
f(R) models in this class. A few examples of the f(R) 
functions are shown in Fig. [TJ 

There is no true cosmological constant introduced in 
this class, unlike in the models of [85j |. However, at cur- 
vatures high compared with m 2 , f(R) may be expanded 
as 



lim f(R) fa -— m 2 



Ci 2 



m 



(7) 



Thus the limiting case of c x jc\ — > at fixed C\ jc^ is a cos- 
mological constant in both cosmological and local tests 
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of gravity, as we shall see. Moreover, at finite Ci/c|, the 
curvature freezes into a fixed value and ceases to decline 
with the matter density, creating a class of models that 
accelerate in a manner similar to ACDM. These models 
therefore also do not exhibit the problems of "mCDTT"' 
models with the form f(R) = p 4 /R. (Note the sign dif- 
ference from the original inverse curvature CDTT model 
While these models can accelerate the expansion, 
they evolve in the future into an unstable regime where 
1 + f R < and also do not contain ACDM as a limiting 
case of the parameter space [82[ . Note that n = 1 resem- 
bles the mCDTT plus a cosmological constant at high 
curvature. Likewise n = 2 resembles the inverse curva- 
ture squared model [HI plus a cosmological constant. 



B. Background Evolution Equations 

Variation of the action ^ with respect to the metric 
yields the modified Einstein equations 



G a p + fRRap — 
where the field, 



/ 



□/A 9afi - = K T af3 ,(8) 



df(R) 
dR 



(9) 



will play a central role in the analyses below. 

Since modifications only appear at low redshift, we 
take a matter-dominated stress-energy tensor. For the 
background Friedmann-Robertson- Walker (FRW) met- 
ric, 



R = 12H 2 + QHH' . 



(10) 



where H(lna) is the Hubble parameter and ' = d/dlna. 
The modified Einstein equations become the modified 
Friedmann equation 



H 2 - f R (HH' + H 2 )+ l -f + H 2 f RR R' 
6 



2 - 
K Z p 



(11) 



To solve these equations, we re-express them in terms 
of parameters the values of which vanish in the high- 
redshift limit where f(R) modifications are negligible 



!JH 
Vr 



H 2 

m? 
R 



— 3a 



-3 



(12) 



Equations (fTT)|) and (TlTj) become a coupled set of ordinary 
differential equations 



1 



Vh 



Vr 



Vr - 4y ff . 



3 

9a~ 3 



Vh 



-3 m 2f 



RR 



(13) 
(14) 



r | 1 1 -3 

Vh - Jr I -Vr -VH-^a 



6 m 2 



To complete this system, we take the initial conditions 
at high redshift to be given by detailed balance of per- 
turbative corrections to R — k 2 p. 

The impact of f(R) on the expansion history can be 
recast as an effective equation of state for a dark energy 
model with the same history 



1 



Weft = — 
3vh 



(15) 



The two equations (fl3|) and (| 14[) combine in the high 
curvature limit — a -3 3> UhiVr — to form 



y'k 



i 



Wh 



3m 2 f RR a 



Z^Uh — driving terms , (16) 



where [. . .] contains a time-dependent friction term whose 
exact nature is not relevant for the qualitative argument 
and the driving terms involve the matter density. As 
shown in [6lj , a critical requirement of the model is 
that f RR > such that Eqn. (fT6|) becomes an oscilla- 
tor equation with real — not imaginary — mass. For the 
background, it is convenient to express this as the di- 
mensionless quantity 



B = 



fRR H_ 



1 + fR H> 



(17) 



This oscillator equation and the parameter B have a sim- 
ple interpretation in terms of the scalar field f R as we 
shall now see. 



C. Field Equations 

The impact of f(R) can be alternatively viewed in 
terms of the field equation for f R . The trace of Eqn. (j8]) 
can be interpreted as the equation of motion for f R 



30f R -R + f R R-2f = -K 2 p. 
This equation can be recast in the form 

dV eS 



□/« = 

with the effective potential 
dV eS 1 



df R 



R-f R R + 2f-K J p) 



(18) 



(19) 



(20) 



df R 3 

The effective potential has an extremum at 

R-Rf R + 2f = K 2 p. (21) 

In the high-curvature regime, where \f R \ -C 1 and 
<C 1, the extremum lies at the general-relativistic 
expectation of R = n 2 p. The curvature at the extremum 
is given by 



in 



fR 



d 2 V cS = I ( 1 + fR 
dfl 3 f RR 



R 



(22) 
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FIG. 2: Cosmological evolution of the scalar field fa and the 
Compton wavelength parameter B for models with n — 1,4. 
Both parameters control observable deviations from general 
relativity and deviations decline rapidly with redshift as n 
increases. 



and hence the extremum is a minimum for B > and 
a maximum for B < in the high-curvature limit with 
\Jr\i \fnRR\ ^ !• Finally the Compton wavelength 



.In 



implies that, in this limit, 



(23) 



(24) 



such that B 1 / 2 is essentially the Compton wavelength of 
fn at the background curvature in units of the horizon 
length. The Compton wavelength plays an important 
role in both cosmological and local tests of f(R) models, 
as we shall see. 



D. Expansion History 

We now evaluate the expansion histories for the class of 
f(R) models in Eqn. ([3]). First, we would like to narrow 
the parameter choices to yield expansion histories that 
are observationally viable, i.e. that deviate from ACDM 
in the effective equation of state (fT5|) by no more than 1 1+ 
w e g\ < 0.2 during the acceleration epoch. This equates 
to choosing a value for the field at the present epoch 
/^o = //i(hia = 0) <SC 1 or, equivalently, Rq 3> m 2 . In 
this case, the approximation of Eqn. (JT]) applies for the 
whole past expansion history and the field is always near 
the minimum of the effective potential 

R = n 2 p - 2/ re n 2 p + 2— m 2 , (25) 

where the 2f term is nearly constant and mimics the 
energy density of a cosmological constant. Thus, to ap- 
proximate the expansion history of ACDM with a cosmo- 
logical constant 57a and matter density Sl m with respect 



to a fiducial critical value, we set 



— ~ 6— A 



C2 



r . 



(26) 



leaving two remaining parameters, n and ci/c 2 , = 
6f2A/c2^m to control how closely the model mimics 
ACDM. Larger n mimics ACDM until later in the expan- 
sion history; smaller ci/c 2 , mimics it more closely. Note 
that, since the critical density and Hubble parameter de- 
pend on the fn modification, Q m is only the true value 
in the limit 



lim n m = Q, m , 

Ci/c|— >o 



(27) 



whereas the matter density in physical units remains un- 
changed Cl m HQ — VL m HQ. 

For the flat ACDM expansion history 




R ps 3m J 



and the field takes on a value of 

-R 



(28) 



(29) 



At the present epoch 



Ro 
Bo 



12 



Cl 



12 



6n(n +1) a f 12 



(i + /Ho)n 



m °2 



f2, 



-9 



-n-2 



(30) 



In particular, for (l m = 0.24 and Ha = 0.76, Rq = 41m 2 , 
fm w -n Cl /c 2 /(41)" +1 and B w -0.61(n + l)f R0 
f° r |/-Ro| *C 1. The consequences of cosmological and 
solar system-tests can be phrased in a nearly model- 
independent way by quoting the field value Jb- Con- 
sequently, we will hereafter parameterize the amplitude 
Ci/c 2 , through the cosmological field value today, //jo- 
in Fig. [2j we show several examples of the background 
evolution of For a fixed present value fno, a larger n 
produces a stronger suppression of the field at high red- 
shift and a larger value of B relative to fa. The steepness 
of this suppression will play an important role for galactic 
tests in €TlTD1 



The effective equations of state for these models are 
shown in Fig. [3] Deviations from a cosmological constant, 
W e ft = —1, are of the same order of magnitude as fuo- 
This class of models has a phantom effective equation of 
state, w e g < — 1, at high redshift and crosses the phan- 
tom divide at a redshift that decreases with increasing 
n. Note that an effective equation of state that evolves 
across the phantom divide is a smoking gun for modified 
gravity acceleration or dark energy with non-canonical 
degrees of freedom or couplings 8t 
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FIG. 3: Evolution of the effective equation of state for n — 1,4 
for several values of the cosmological field amplitude today, 
fuo- The effective equation of state crosses the phantom di- 
vide W e a = — 1 at a redshift that decreases with increasing n 
leading potentially to a relatively unique observational signa- 
ture of these models. 



E. Linear Perturbations 
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FIG. 4: Fractional change in the matter power spectrum P(k) 
relative to ACDM for a series of the cosmological field ampli- 
tude today, /ho, for n = 1, 4 models. For scales that are below 
the cosmological Compton wavelength during the acceleration 
epoch k > (aH)B 1 ' 2 perturbation dynamics transition to the 
low-curvature regime where 7 = 1/2 and density growth is 
enhanced. This transition occurs in the linear regime out to 



field amplitudes of \fh 



- 10" 



Given an expansion history that defines fn(lna) and 
S(lna), the evolution of linear perturbations can be 
solved using the techniques of [6l|. The principal fea- 
ture of the linear evolution is that once the wavelength 
of the perturbation becomes smaller than the Compton 
wavelength in the background 



aH 



-.B 1 / 2 > 1 



(31) 



strong deviations from the GR growth rate appear. In 
particular, the space-space $ and time-time 'J pieces of 
the metric fluctuations in the Newtonian (longitudinal) 
gauge evolve to a ratio 



1 



* 7 2 ' 



(32) 



implying the presence of order- unity deviations from GR. 

The consequence of this relative enhancement of the 
gravitational potential ^ is an increase in the growth rate 
of linear density perturbations on scales below the Comp- 
ton wavelength. If the Compton wavelength is longer 



than the non-linear scale of a few Mpc, this transition 
leads to a strong and potentially observable deviation in 
the matter power spectrum. Even percent-level devia- 
tions in the power spectrum are in principle detectable 
with future weak-lensing surveys. If the Compton wave- 
length approaches the horizon, it can substantially alter 
the CMB power spectrum as well [6l|. 

In Fig. [U we illustrate this effect for n = 1 and n = 4 
models. Deviations occur in the linear regime down to a 
field amplitude of \fu\ ~ 10~'. For these small field 
amplitudes, the expansion history and hence distance 
measures of the acceleration are indistinguishable from 
a cosmological constant with any conceivable observa- 
tional probe. Nonetheless, linear structure can provide a 
precision test of gravity that, we shall see, rivals that of 
local tests in a substantially different curvature regime. 
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III. LOCAL TESTS 

In this section, we consider local tests of f(R) gravity. 
We begin in mil Al with a general metric around spher- 
ically symmetric sources and its relationship to the f R 
field. In §111 B| we discuss the qualitative behavior of 
the field solutions and their relationship with the Comp- 
ton wavelength. We evaluate solar-system constraints in 
§111 CI and the requirements they place on the extent and 
evolution of the galactic halo in §111 Dl 



A. Metric and Field Equations 

We take the general spherically symmetric isotropic 
form for the metric around a source centered at r = 

ds 2 = -[l-2A(r) + 2B(r)]dt 2 

+ [l + 2A(r)](dr 2 +r 2 dn), (33) 

where we assume that \A(r)\ -C 1 and \B{r)\ -C 1 near 
the source, such that the metric is nearly Minkowski. In 
the GR limit, B(r) — > 0; limits on B in the solar system 
provide the strongest tests of modification to gravity of 
the type considered here. 

By definition of the Ricci tensor, the sources of the 
metric potentials A and B are given by 



V 2 (A + B) = -±R, 

V 2 B = -\{R\ + R/2) 



(34) 
(35) 



As an aside, choosing models for which B > in mi Al 
not only results in the existence in the expansion history 
of a stable matter-dominated era, but ensures that the 
models do not exhibit the related instability for stellar- 
type objects [§(| HH- Small, time-dependent perturba- 
tions to a high-curvature solution of Eqn. (|39p have pos- 
itive mass squared and do not grow in this class of f(R). 

In the limit that \f R \ < 1 and \f/R\ <C 1, valid for all 
sources that we shall consider, 



V 2 A « -~k 2 P +^(k 2 P -R) 
V 2 B « i( K 2 P -i?) . 



(40) 
(41) 



As expected, the source of B is the deviation of the cur- 
vature R from the GR value of n 2 p. Moreover, in the 
same limit, equation (|39p for f R becomes 



V 2 /*« -{R-n 2 p) 



(42) 



A solution for f R therefore gives B up to constants of 
integration 



B(r) = -f R (r) + ai + 



a 2 



(43) 



Since B must remain finite at r — 0, <22 = 0. Let us 
assume that at sufficiently large radii f R (r) — » f Roc and 
B then 



B(r) = -\f R {r) - fn, 



-A/*(r) 



(44) 



Note that a low-curvature R <§; K 2 p solution may be 
accommodated byB~ —A. 

To relate the time-time component of the Ricci tensor 
to f(R), let us take the time-time component of the field 
equation JHJ 

(1 + f R )R° -l(R + f) + nf R + dff R = ~n 2 p . (36) 



For radii beyond which the source K 2 p — R becomes neg- 
ligible B(r) cx 1/r. It is convenient to then define an 
effective enclosed mass 



such that 



Af cff = 4tt (p- R/n 2 )r 2 dr , 



(45) 



Under the assumption of a static solution, we can com- 
bine this equation with the trace equation (|18[) to obtain 



R n = 



-2n 2 p+\R-\f + f R R 
3(1 + f R ) 



(37) 



Even in a low-curvature solution where R -C K 2 p, R° — 
0(n 2 p). 

Eqn. (|35|) then becomes 



V Z B = -- 



1 / —An 2 p + AR + 5f R R — f 



3(1 + In) 



, (38) 



where f(R) is given by the solution to the trace equation 
(fT8l) in the static limit 



3V 2 f R -R + f R R-2f 



-K 2 p . 



(39) 



B(r) = -Af R (r) 



2GM cS 
3r 



(46) 



Note that this is an implicit solution, since R(f R ). Nev- 
ertheless, we shall see in the next section that it sheds 
light on the behavior of explicit solutions. 

Finally, the deviation from the general-relativistic met- 
ric is given by 



7-1 = 



B 



2M, 



cir 



A-B 



3M t , 



M, 



err 



(47) 



where M tot is the total mass of the system. The two 
limiting cases are M c g <C M to t for which 7 — 1 = 
-2M eff /3M tot and M cS = M tot for which 7 - 1 = -1/2. 
These solutions correspond to high curvature, R ~ K 2 p, 
and low curvature, R <C n 2 p [zl, respectively. 
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B. Compton and Thin-Shell Conditions 

Before examining explicit solutions for /r(t") given 
p{r), we show how the nature of the solutions is tied 
to the Compton wavelength of the field and exhibits the 
so-called chameleon mechanism [78[ for hiding scalar de- 
grees of freedom in high-density regions [79|, [8(| • 

An examination of Eqn. ([39]) shows that there are 
two types of local solutions to the field equations dis- 
tinguished by a comparison of the Compton wavelength 
of the field fa 



RR ■ 



(48) 



to the density structure of the source. 

Let us again assume that <C 1 and \ f/R\ <§C 1 so 
that the field equation is well approximated by Eqn. (|4"2")l . 
The first class of solutions to this equation has high cur- 
vature R w n 2 p and small field gradients V 2 /^ <C n 2 p. 
The second class of solutions has low curvature i? <§; n 2 p 
and large field gradients V 2 /^ « —n 2 p/3. 

A sufficient condition for the high-curvature solution 
is that field gradients can be ignored at all radii when 
compared with the density source. A necessary or con- 
sistency condition is that field gradients implied by the 
high-curvature solution f(R = n 2 p) can be ignored com- 
pared with local density gradients. More specifically 



Jrr 



R=K 2 P 



d 2 p < p , 



,1/2 

■Irr 



R=k 2 p 



d iP <p, (49) 



i.e. that the density changes on scales that are much 
longer than the Compton wavelength. Equivalently, a 
mass source induces changes in the field with a Yukawa 
profile of e~ mf n r /r which are highly suppressed on scales 
larger than the Compton wavelength. We call this con- 
dition the Compton condition. 

If the Compton condition is satisfied at all radii, then 
the high-curvature solution is also valid at all radii and 
deviations from GR will be highly suppressed. If this con- 
dition is violated beyond some outer radius, then a por- 
tion of the exterior must be at low curvature R <C n 2 p. 
Moreover, since Birkhoff's theorem does not apply, the 
exterior low-curvature solution can penetrate into the re- 
gion where the Compton condition ([4"9"| is locally satis- 
fied. The interior solution then depends on the exterior 
conditions. We have seen in §11 El that linear cosmological 
perturbations are in the low-curvature regime on scales 
smaller than the cosmological Compton wavelength and 
so the Compton criteria must be violated far in the ex- 
terior if \/ro\ > 10~ 7 . We will return to this point in 

To quantify these considerations, note that the max- 
imal change in fa from the interior to the exterior is 
imposed by the low-curvature assumption R <C n 2 p or 

M eff = M tot 



A/fl(r) < ^mW, 



(50) 



where $Af(^) is the Newtonian potential profile of the 
source, i.e. $m ~ GM tot /r exterior to the dominant 
mass. This condition sets an upper limit on the difference 
between the interior and exterior field values for a static 
solution. 

If the thin-shell condition is satisfied and \Afn(r)\ <C 
^Ai(r), then M 3> M e ff and somewhere in the interior 
there must exist a high-curvature region where R — > k 2 p. 
To estimate where this occurs consider a local version of 
Eqn. (J50J) 



^ « 2 (p- 



Poo)r 2 



(51) 



From the outside in, when this condition is first satisfied, 
there is enough source to make the transition between the 
interior and exterior field values. Once this is satisfied, it 
remains satisfied in the interior as long as further changes 
in fji are much smaller than the initial jump. In other 
words, the exterior field is only generated by the "thin 
shell" of mass M c g that lies outside of this transition. We 
will call this the thin-shell criterion and such a solution 
is known in the literature as a chameleon solution. 

The thin-shell criterion is related to the low-curvature 
linearization condition of [77 ] Eqn. (13). There, a solu- 
tion for the curvature is found by linearizing Eqn. (|39]) 
around its background value. Requiring that the value of 
the perturbation be smaller than that of the background 
curvature, results in exactly the opposite of Eqn. (|50p 



|/i?.o 



> 3<Mr). 



(52) 



Therefore, the linearization procedure is not valid for ex- 
actly those sets of parameters for which the thin-shell 
condition is satisfied: high-curvature solutions are nec- 
essarily non-linear. When the linearization is valid, the 
solution is low-curvature everywhere leading to large de- 
viations from GR in the interior. 

The thin-shell criterion (jSTj) is also related to but 
stronger than the local Compton condition 



Irr < r p , 



(53) 



where r p is the distance over which the density field 
changes. Converting derivatives to finite differences 
Irr « Sf R /SR 



SfR < SRr 2 p « K 2 5pr 2 p . 



(54) 



The difference between the two conditions is that the 
Compton condition involves the small change in the field 
<5/t?. a t high curvature, whereas the thin shell criterion in- 
volves the potentially larger change in the field Afn from 
the high- to low-curvature regimes. If there is no transi- 
tion to low curvature in the exterior then these conditions 
are the same. 

The thin-shell condition implies that the field does not 
always sit at the local potential minimum R w k 2 p. 
Nonetheless, the field does choose the energetically fa- 
vorable configuration. The field will not lie at the 
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FIG. 5: Density profile in the solar interior and vicinity (solid 
curve). Under general relativity (GR), the curvature R would 
track the density profile. For the f(R) model with n = 4, 
cosmological field amplitude |/no| =0.1 and a galactic field 
amplitude that minimizes the scalar potential, the curvature 
tracks the GR or high-curvature limit out to the edge of the 
solar corona at about 1AU (dashed line). 

potential minimum in the interior region if the energy 
cost for introducing a field gradient between the interior 
and exterior is too high. The potential energy density 
cost for not lying at the potential minimum is 

A*/ „ ^ p \Af R \ , (55) 

where A/r is the difference between potential minimum 
and the exterior solution. Compare this with the gradient 
energy density gain from not introducing a profile in the 
field 

Ai?~I^M. (56) 
2 r l 

The potential energy cost outweighs the gain if 

|A/fl| < n 2 pr 2 , (57) 
which is the thin-shell condition. 



C. Solar System to Galaxy 

We now consider explicit solutions of the field equa- 
tion (|39[) . given the density profile p(r), in the solar vicin- 
ity. In this section, we will assume that the galaxy has 
sufficient mass and extent to bring the field to the poten- 
tial minimum of Eqn. ([29]) in the outskirts of the solar 
system. In the next section, we will discuss the require- 
ment this assumption places on the structure and evolu- 
tion of the galactic halo. 

Specifically, we set the boundary condition /.Roo = fug 
where 

fRg = MR = K 2 Pg ) . (58) 



FIG. 6: Solar solution for fn(r) with n = 4 for a series of 
cosmological field amplitudes fao with a galactic field that 
minimizes the potential. The field smoothly transitions from 
the fa ~ interior value to the galactic values once the thin- 
shell criterion is violated. 



Here, p g is the average galactic density in the solar vicin- 
ity. Under this assumption, the galactic field value is 
related to the cosmological one as 

\f m J U m 2 0.13 Vl0- 24 gcm-3j 

(59) 

For the density profile in the solar vicinity we take 
the helioseismological model of [g^] for the solar interior, 
the chromosphere model of (93| and the corona model of 
[qH for the surrounding region. We add to these density 
profiles a constant galactic density of p g = 1CU 24 g cm~ 3 . 
This profile is shown in Fig. O (solid curve). 

We solve the field equation |39|) as a boundary-value 
problem using a relaxation algorithm. The source-free 
(p = 0) field equation has exponentially growing and 
decaying Yukawa homogeneous solutions exjp(±mf R r)/r. 
Initial- value integrators have numerical errors that would 
stimulate the positive exponential, whereas relaxation 
methods avoid this problem by enforcing the outer 
boundary at every step. We show an example solution in 
Fig. [H (dashed curve) for n = 4 and \fm\ =0.1. 

A solution is found on an interval, by requiring that the 
field fn minimize the potential (|20[) at two chosen radii: 
one far away from the Sun and its corona, the other — 
inside the solar density distribution. We place the outer 
boundary at r — lO 6 r0, where the density distribution is 
entirely dominated by the constant galactic-density com- 
ponent. The solutions are robust to increasing the radius 
of this boundary. 

For the inner boundary, we take the starting point as 
approximately 1000 Compton wavelengths from the tran- 
sition to low curvature; the Compton condition is well 
satisfied there. Interior to this point, the solution is more 
efficiently obtained by a perturbative solution around the 
high-curvature solution of RS ' — n 2 p. The first-order 
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correction to this solution is R — RS°> + R^ with 
R {1) = (3W 2 f R + f R R-2f) 



(60) 



With this correction to R, the first-order correction to 
fa can be obtained using initial-value-problem meth- 
ods. This series can be iterated to arbitrary order and is 
strongly convergent when the Compton condition is sat- 
isfied. In fact, at our chosen starting point, the zeroth or- 
der solution typically suffices. We have also checked that 
for cases where the Compton condition is everywhere sat- 
isfied, the numerical solution matches the perturbative 
solution. Finally, we stop the relaxation code when the 
solution satisfies the trace equation |39|) to at least ICP 6 
accuracy. 

Let us now relate the numerical solutions shown in 
Fig. [6] to the qualitative analysis of the previous section. 
First, consider the Compton condition. The Compton 
wavelength at R 



n 2 p is 



\ fR = (10.6pc) (8.14 x 10 

x[(" + I)|/™|] 1/2 
P_ 

-24 



-7\ («-l)/ 2 



R tt m h 2 
ri* 0.13 

(n+2)/2 



2\ (n+l)/2 



I0- 24 g cm" 

For example, for n = 4 and the fiducial cosmology 

-3 



X fR « (83OOr )|/ flo | 1/2 



10" 24 g cm" 3 



(61) 



(62) 



so that the Compton condition is satisfied for the whole 
solar profile for \fm\ < 10~ 2 . Note that, even at the base 
of the corona where the density is p w I0~ 15 g cm -3 , 
the Compton wavelength is ~ lO _23 |/flo| 1//2 r0. Thus, 
in spite of the steep density gradient at the edge of the 
Sun, the Compton condition is well satisfied, allowing the 
solution in Fig. [5] to follow the high-curvature solution 
R = K 2 p closely through the transition. 

On the other hand, the thin-shell criterion is satisfied 
in the solar corona up to |/,ro| < 10 _1 . Thus for n = 4, 
we expect order-unity cosmological fields to be achiev- 
able with the entire solar interior including the edge and 
chromosphere in the high-curvature regime. 

The numerical solutions shown in Fig. [S] verify the 
qualitative behavior described in the previous section. 
For |/ijo| < 10 -2 , the deviations from the high-curvature 
R = K 2 p limit are fractionally small since the Compton 
condition is everywhere satisfied. For |/,ro| < lO -1 , the 
break to low curvature occurs in the corona. This break 
occurs gradually in the field profile fn{r) but rapidly in 
the curvature (see Fig. [7|. At small field values, a small 
change in f R represents a large change in the curvature. 
M e s is approximately just the mass between this transi- 
tion and the point at which the galactic density exceeds 
the corona. Outside of this transition, the exterior field 
relaxes to the galactic value as |A/r| oc e~ mf R r jr. In 
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FIG. 7: Solar solution for the curvature R for the n = 4 
model and a series of cosmological field amplitudes Jrq with 
a galactic field that minimizes the potential. The solution 
abruptly transitions from the high-curvature R ~ K 2 p to the 
low-curvature regime once the thin-shell criterion is violated. 



these examples, the Compton wavelength in the galaxy 
is of order 10 3 — 10 4 r Q and the mass term further sup- 
presses the deviations from GR. 

In Fig. [8] we show |7 — 1| for the same n — 4 models. 
The deviations peak at ~ I0~ 15 . Such deviations easily 
pass the stringent solar system tests of gravity from the 
Cassini mission 19 511 



|7- 1| < 2.3 x 10" 



(63) 



under the assumption that the galactic field fn g is given 
by the potential minimum. 

Models that saturate this observational bound have a 
sufficiently large that the thin-shell criterion is first 
satisfied at the edge of the Sun, where the fractional en- 
closed mass becomes 2M c s/3M to t ~ 10 -5 (see Eqn. |47|) . 
For a given f R o, this can be achieved by lowering n. In 
such models, the additional contribution to the effective 
mass outside of the photosphere is negligible and the field 
solution obeys Eqn. (|46|) . The exterior solution therefore 
becomes 



A/*(r) w (7 - l)*w(r) w (7 - 1) 



GM P) 



(64) 



since the Compton wavelength in the exterior implied by 
Eqn. (|61[) is much larger than the solar system. 
Given that 



GM, 



r@ 



2.12 xIO" 6 . 



(65) 



the solar-system constraints can be simply stated as 

|A/*(r )|«|/^|< 4.9xl0- n . (66) 

Note that this bound is independent of the form of 
f(R) and the assumption that fn g is given by the min- 
imum of the effective potential. The model dependence 



10 



1.5 



o 



0.5 - 





i-i—A - 
ft — L r 


lf„J-0 0*1 




- — !fW=o.oi 




1^,1=0.001 




; — Lf fi0 i=o / 








y 







100 



1000 

rlr a 



io 4 



FIG. 8: Metric deviation parameter I7 — 1| for n — 4 mod- 
els and a series of cosmological field amplitudes fm with a 
galactic field that minimizes the potential. These deviations 
are unobservably small for the whole range of amplitudes. 



comes from the implications for the cosmological field 
value fno. Using Eqn. (|59p at a galactic density of 
p g = 10~ 24 g cm -3 , we can translate this into a bound 
on the amplitude of the cosmological field 



I/hoI < 74 (1.23 x 10 6 ; 



71-1 



Rq fl m h 2 
mP 0.13 



-(ti+1) 



(67) 



As shown in Fig. [5] this is a fairly weak constraint. For 
11 > 1 it allows order unity cosmological deviations from 
GR. Note that the models of [8(| are equivalent to a 
continuation of the approximate form of our model in 
Eqn. but with n < 0. Solar system constraints on 
cosmological amplitudes are significantly stronger in that 
class. Finally, although a detailed calculation is beyond 
the scope of this work, laboratory constraints on fifth 
forces are also weak under the same assumptions given 
the much larger effective densities involved [781 ] . 



D. Galaxy to Cosmology 



The cosmological bound Eqn. ([67)) from the solar- 
system constraint Eqn. (|66| is robust but weak. A related 
but potentially more powerful constraint comes from the 
transition between the high curvature of the galaxy and 
the desired cosmological curvature. Indeed, the implicit 
assumption in applying the solar-system constraint to the 
cosmology is that the galaxy itself is in the high-curvature 
regime with respect to its own density profile. 

The validity of this assumption depends on both the 
structure of the galactic halo and its evolution during 
the acceleration epoch. Large |/ro| requires that ei- 
ther the galactic gravitational potential is substantially 
deeper than in ACDM or that it has not yet reached its 
equilibrium value. 

From the linear theory analysis in ijllEl we know that 
the Compton condition is violated for linear perturba- 




FIG. 9: Maximum cosmological field \fm\ allowed by the so- 
lar system constraint alone (top curve) under the assumption 
that the galactic field remains at the potential minimum at 
the present epoch. Models with \fm\ > 10 -6 (bottom curve) 
have galactic fields |/a 9 | that evolve to higher values dur- 
ing the acceleration epoch, potentially enabling substantially 
stronger constraints with structure formation simulations. 



tions if |/_ro| > 10 -7 . Hence, at the outskirts of the galac- 
tic halo where the density profile joins onto the large-scale 
structure of the universe, there must be a transition from 
high to low curvature. As we have seen in ^IIIBl the low- 
curvature cosmological field will eventually penetrate into 
the galaxy unless 



\Af R (r)\ = \f Rg -f m \tt\f m \<-$ g 



(68) 



where <& 9 is the Newtonian potential of the galaxy. For 
definiteness, let us consider the NFW density profile 



Pg( r ) 



1 



47r r(r 



(69) 



where M g is the galactic mass contained within 5.3r s 
and r s is the scale radius of the dark-matter halo. This 
density profile has a Newtonian potential of 



GM a 



ln(l + r/r s ) 



and a maximum rotation-curve velocity of 



0.46 



(70) 



(71) 



at 2.16r s . Taking the thin-shell criterion to be satisfied 
at r sa r s , such that the interior is in the high-curvature 
regime, leads to an upper limit on |/.ro| for a static solu- 
tion of 



l/reo) <2x 10 



-°(- 



V300km/s 



(72) 



For higher cosmological values of the field |/ro|, the 
galaxy will relax over time to the low-curvature solu- 
tion to minimize the cost of field gradients. Note that 
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FIG. 10: A galactic solution for the scalar curvature R for 
n = 1/2 with an NFW density profile (|69jl and conserva- 
tive parameters r s — 75 kpc and v maK = 300 km s _1 . For 
I /ho | ^ 2 x 10~ 6 , no static solution can be found which has 
high curvature, R » K 2 p, inside the halo. As in the solar case, 
once the thin-shell criterion is violated the solution abruptly 
transitions to low curvature and matches onto the cosmolog- 
ical value of /rq outside the halo. 



even a value of \fm\ = 10~ 6 provides potentially observ- 
able modifications in the linear regime since gravity at 
the juncture is still modified by order unity (see Fig. HJ). 
We have verified through the numerical techniques of the 
previous section that the constraint in Eqn. (|72[) is nearly 
independent of the functional form of f(R). In Fig. ITQ1 
we show an example solution. Here we err on the conser- 
vative side by taking v max — 300 km s _1 and r s — 75kpc 
to reflect a somewhat more massive and extended halo 
by a factor of a few than expected in ACDM for our 
galaxy (c.f. (96|, Fig. 7). Empirical data exists only 
out to ~ 20kpc where the rotational velocity reaches 
v « 230 km s" 1 . 

Unfortunately, the bound in Eqn. (|72")l is suggestive 
but not definitive. A cosmological simulation will be re- 
quired to determine how haloes of galactic size, which 
are embedded in a group-sized dark-matter halo, which 
itself is part of the quasi-linear large-scale structure of 
the universe, evolve during the acceleration epoch. The 
density profiles of the structures in which the galaxy is 
embedded can further shield the galactic interior from 
the low-curvature solution. Furthermore, the cosmologi- 
cal background itself was at high curvature at z ^ 1 . For 
n ^> 1, the condition that the background \fu\ < 10~ 6 is 
well satisfied for z » 1 even for approaching unity. 
For example if n — 4, the cosmological field drops by 10 3 
by z = 2 (see Fig. [2]) . Therefore when the galactic halo 
formed both its interior and exterior were at high cur- 
vature. The local curvature would then follow the local 
density R rs k 2 p closely everywhere. 

Only during the recent acceleration epoch is the thin- 
shell criterion for the galaxy violated at \fm\ > 10~ 6 . 
The low curvature, high field values of the background 



will then begin to propagate into the interior of the galac- 
tic halo in a manner that requires a simulation of the 
process. The static, thin-shell bound for the galaxy given 
by Eqn. (f72"|) is therefore overly restrictive but suggests 
that cosmological simulations should enable much more 
stringent bounds on f(R) models than solar-system test 
alone. 



IV. DISCUSSION 

We have introduced a class of f(R) models that ac- 
celerate the expansion without a cosmological constant. 
Its parameters allow the gravitational phenomenology ex- 
hibited in cosmological, galactic and solar-system tests to 
span the range between infinitesimal and order-unity de- 
viations from general relativity with a cosmological con- 
stant. 

In these models, unlike the original model of [l[ and 
related generalizations, the general-relativistic or high- 
curvature value for the Ricci scalar R w K 2 p is the so- 
lution that minimizes the potential for the scalar degree 
of freedom fn = df(R)/dR. This feature is critical for 
both cosmological and solar-system tests. 

Solar-system tests of f(R) gravity alone place only 
weak bounds on these models, despite a strong and nearly 
model-independent limit on the /r field amplitude in the 
galaxy of < 5 x 10~ n . Indeed, we show that for 

a range of models that include cosmological fields of or- 
der unity \fjta\ < 1, the field deviates from the high- 
curvature regime only for a brief interval in the solar 
corona. Likewise, deviations from general relativity in 
the metric are generated by the mass in this interval 
and so are unobservably small. This is strikingly differ- 
ent from constraints that would arise if one embeds the 
Sun directly into a medium at the cosmological density 
\f R0 \ < 5 x 10- 11 . 

Without further constraints on the size and evolution 
of the galactic halo, solar-system tests may be evaded 
relatively easily. This is because of the assumed strong 
density scaling of the field amplitude at the potential 
minimum. However, though the galactic field begins at 
high redshift at its potential minimum, it will only remain 
there if the galaxy is sufficiently massive to protect it 
against the cosmological exterior which is evolving to low 
curvature at z < 1. 

An order of magnitude estimate, based on the extrap- 
olation of rotation-curve measurements and the assump- 
tion that the galactic halo does not differ substantially 
from ACDM expectations in the outskirts, suggests that 
an isolated galaxy that is otherwise like our own will only 
remain stably at high curvature if the cosmological field 
18 l/flol 10 -6 . Despite the fact that the high-curvature 
solution still minimizes the potential energy of the field, 
the gradient energy implied by the field profile from the 
galactic interior to exterior is too high. This estimate is 
also nearly independent of the functional form of f(R) 
and improvements in the solar-system constraint. Turn- 
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ing this estimate into a firm constraint on models will 
require cosmological simulations of f(R) acceleration to 
examine how far into the solar system the exterior low- 
curvature field values can penetrate by the present epoch 
in the local environment of our galactic halo. 

Distance-based measurements of the expansion history 
will be limited to testing |/ fl0 | > 10~ 3 — 10~ 2 for the 
foreseeable future, since the field amplitude determines 
the deviations in the effective equation of state to be 
of comparable size. Nonetheless, future, percent-level 
constraints on the matter power spectrum in the linear 
regime offer potentially even stronger tests of f(R) mod- 
els than solar or galactic constraints, in principle down 
to amplitudes of \fm\ ~ 1CP 7 . This sensitivity is due to 
the large Compton scale in the background across which 
perturbations make the transition from low to high cur- 
vature and exhibit order-unity deviations from general 



relativity. Cosmological simulations are also required to 
determine how these signatures can be disentangled from 
the non-linear evolution of structure. 
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